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ABSTRACT 

The a-formalism is a common way to parametrize the common envelope interac- 
tion between a giant star and a more compact companion. The a parameter describes 
the fraction of orbital energy released by the companion that is available to eject the 
giant star's envelope. By using new, detailed stellar evolutionary calculations we de- 
rive a user-friendly prescription for the A parameter and an improved approximation 
for the envelope binding energy, thus revising the a equation. We then determine a 
both from simulations and observations in a self consistent manner. By using our own 
stellar structure models as well as population considerations to reconstruct the pri- 
mary's parameters at the time of the common envelope interaction, we gain a deeper 
understanding of the uncertainties. We find that systems with very low values of q (the 
ratio of the companion's mass to the mass of the primary at the time of the common 
envelope interaction) have higher values of a. A fit to the data suggests that lower 
mass companions are left at comparable or larger orbital separations to more massive 
companions. We conjecture that lower mass companions take longer than a stellar 
dynamical time to spiral in to the giant's core, and that this is key to allowing the 
giant to use its own thermal energy to help unbind its envelope. As a result, although 
systems with light companions might not have enough orbital energy to unbind the 
common envelope, they might stimulate a stellar reaction that results in the common 
envelope ejection. 

Key words: binaries: close — planetary nebulae: general — stars: horizontal branch 
— stars: evolution — stars: white dwarfs 



1 INTRODUCTION 

Common envelope (CE) binary interactions occur when ex- 
panding stars transfer mass to a companion at a rate so high 
that the companion cannot accrete it. This results in the 
companion being engulfed by the envelope of the primary 
jPaczvnskillllffli ). The companion's orbital energy and an- 
gular momentum are then transferred to the envelope via 
an as yet poorly characterised mechanism. This can result 
in the ejection of the envelope and in a much reduced orbital 
separation. If the companion cannot eject the CE, it merges 
with the core of the primary. The CE interaction is thought 
to last for only a few years, and it is therefore likely that we 
have never witnessed it dire ctly, although some cl aims have 
been made, e.g., for V Hya (|Kahane fc Jura 199f|). The ex- 
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istence of companions in close orbits around evolved stars, 
whose precursor's radius was larger than today's orbital sep- 
aration, vouches for such interactions having taken place. 

The CE interaction is thus thought to be responsi- 
ble f or short period binaries such as cataclysmic variables 
(CV; iKindl 19881 : IWarnerlll995t). close bi nary central stars of 
plan etary nebula (PN: |Pe Marcdl2009h. subdwarf B bina- 
ries l|Morales-Rueda et al.ll2003l: iHan etal1l2002l . 120031 ). low 
mass X-rav binaries (|Charles fc CodfeOQq), the progenitors 
of Type la supernovae ( Bclczvns ki et all 120051 ) and other 
classes of binaries and single stars t hought to have suffere d 
a merger (such as FK Comae stars; iBopp fc Stencellll98ll ). 
The specific characteristics exhibited by these binary classes, 
as well as their relative population sizes are dictated by the 
period and mass ratio distribution of the progenitor binary 
population, as well as by the details of the physical interac- 
tion during the CE phase. 
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The CE interaction can be parametrised in te rms of the 
binding and orbital energy sources at play (e.g., IWebbinkl 
1 19841 . l2008l ). and the post-CE period has been expressed 
as a function of how efficiently orbital energy can be used 
to unbind the CE. The efficiency parameter, a, was thus 
introduced: 



-Ebin 
AiJorh 



(1) 



where E^n is the gravitational binding energy of the en- 
velope and AE or ), is the amount of orbital energy released 
during the companion's in-spiral. The expressions used for 
binding and orbital energies in the literature have varied. 
As a result, the conclusions reached in the numerous papers 
discussing the CE interaction by means of the a-formalism 
are difficult to compare. The first motivation of this paper 
is therefore to choose a formalism b y revisiting pas t c hoices . 

Several papers (e.g.. iMaxted et all 120061; 
lAfsar fc Ibanoglul 120081 ; IZorotovic et al.l l201Ch use in- 
dividual post-CE binaries to derive a. Their observed 
primary and secondary masses, together with their orbital 
periods, provide us with parameters of the post-CE systems. 
Based on the primary mass, one can reconstruct the mass 
and radius of the primary at the time of the CE interaction 
and, with this information, a value of a can be derived. 
However, this method has many hidden uncertainties, and 
the values of a derived in this way are, once again, only 
indicative. We therefore use our preferred a-formalism to 
re-evaluate the value of a for observations and simulations 
in a homogeneous way and with an increased attention to 
the sources of uncertainty. 

The values of a determined by simulations, (e.g., 
ISandquist et al.l 1 19981 ) have their own flavour of hidden 
caveats and are not easily comparable with those derived 
from observations. We therefore use the better-understood 
simulations from the literature to gain insight in how the 
values of a determined in this way compare with those de- 
rived from observations. 

In §[5] we discuss the a-formalism in the literature and 
derive our preferred form. We also discuss the value of A, 
often used in parametrizing the envelope binding energy. In 
§ we calculate the value of a for a set of simulations and 
observations in a self consistent manner. We then (§ [4]) dis- 
cuss the dependence of a on stellar and system parameters. 
We conclude and summarise in §[S] 



2 THE a EQUATION 

In this section we determine the best form of the a equation 
(Eq.EQl. 



2.1 The a-formalism in the literature 

The original a equat ion can be found 
iTutukov fc Yungelsonl (1 19791 ). who used: 



, Ir M 2 = M Hc m 
R 2i?Hc ' 



(2) 



where we have maintained the original symbols to empha- 
sise the subtly different assumptions made in each equation. 
/3 is the binding energy parameter equivalent to a, Mae and 



Rue are the post-CE primary's core mass and radius, re- 
spectively, m is the companion's mass, M and R are the 
mass and radius of the primary at the time of the CE in- 
teraction. The value of ^ was taken to be 0.5 to account 
for the fact that the radius of the primary at the start of 
the dynamically- significant part of the CE interaction was 
thought to be twice as large as it was at the beginning of 

the CE in teraction. 

Later llben fc Tutukovl (1 19841 ) used a similar expression 
but modified the symbols: 



Mi 
~A~ 



M 1M M 2 



A 



f 



(3) 



where Mi is the primary mass at the time of the CE and 
Mi,h is the primary (remnant) mass after the CE interac- 
tion. Remembering that there is a factor of 1/2 on both sides 
of Eq.|3]we see that A, the initial binary separation, has re- 
placed the initial primary radius (a reasonable assumption) 
and that Af, the final binary separation, has taken the place 
of the primary core radius, a choice that seems sensible, as 
the final binary separation should be larger than the pri- 
mary's core radius. 

In the same year, IWebbinkl l| 19841 ) rewrote the expres- 
sion as: 



G 



Mi Mi 

\Rl,L 



-acsG 



M lc M 2 MiM 2 



(4) 



2Af 2A t 

where we have adopted the symbols used by IWebbinkl 
(|2008l ): Mi, c and Mi >e are the primary's core and envelope 
masses, respectively, A% is the inital binary separation, Ri : l 
is the Roche lobe of the primary at the onset of mass trans- 
fer and where A is a number of order unity which depends on 
the mass distribution of the primary's envelope (AR is effec- 
tively the mass- weighted mean radius of the envelope) . The 
symbol A was actua lly introduced b y Ide Kool et alj (|l987l ) 
as an addition to the Wcbbink ( 1984) equation and has been 
included in the equation ever since. Eq. [4] contains a term 
for the orbital energy at the beginning of the CE interac- 
tion, a term that others have neglected, on the ground that 
its value is far smaller than that of the final orbital energy, 
due to the considerable in-spiral that happens during the 

inte raction. 

lYungelson et al.l (|l993l ) rewrote the a equation in the 
following way: 



(Mi - Mi, R )(Mi + M 2 ) 
A 



= -aMi R M 2 



(---) 
Af A } 



, (5) 



where all symbols have been previously defined and Ao is 
the binary separation at the beginning of the dynamically- 
significant part of the CE interaction. Eq. [5] actually con- 
tains a factor of 0.5 on both sides of the equality, which was 
cancelled out by lYungelson et al.l (|l993l ). but that must be 
there for the expression of the orbital energy to be correct. 
This factor implies that the size of the primary at the start 
of the dynamically-significant part of the CE interaction is 
2 x Ao (Eq.[3]). Using this expression wil l therefore r esult in 
smaller values of a (as also noticed by Irian et aL I l|l995h ). 
Aside from the d ifference in prim ary radius, Eq. [5] is differ- 
ent from that of IWebbinkl l|l984h also because the binding 
energy of the envelope is more negative: the giant's envelope 
during the CE is bound also by the gravitational attraction 
of the companion within it. 
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Preemptying our derivation in § 12,21 in this work we 
will use: 



G 



XR 



= -aG 



M C M 2 (M c + M e )M 2 



2 At 



2Ai 



(6) 



where M c and M e are now the giant primary's core mass 
(which is assumed to be the same as the mass of the pri- 
mary after the CE interaction) and the primary's envelope 
mass (which is assumed to be ejected by the interaction), 
respectively. In § 12.31 we will use the virial theorem to show 
that the energy budget should include a factor of 0.5 on 
the left-hand-side of Eq. [6] corresponding to the thermal 
energy of the gas lessening the stellar envelope's gravita- 
tional potential well. For now, however, we have presented 
the a equation using only the traditionally-included energy 
sources. The thermal energy source will be discussed further 
in §[231 §E3]and §0 

In the case of observed systems (§ 13. 2[) . we assume 
that the primary is filling its Roche lobe radius at the be- 
ginning of the CE interac tion (i.e., R = Ri.rl) and th at 
At = Ri,RL/n, R L, where l|Egrietor]|l983l ; IWebbinldl200ct l: 



0.49g 

ri ' RL ~ 0.6g~ 2 / 3 + ln(l + g-V3) { '> 

a nd q = Mq/Mfl Our expression is almost identical to that 
of IWebbinkl |l98J), except that our binding energy is some- 
what lower. 

Clearly, irrespective of which formalism one uses, it is 
paramount that values of a calculated with one formalism 
not be used as input to a different formalism in order to 
derive, for instance, values for the final orbital separation. 
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2.2 The binding energy term 

Several authors (e.g.. [Han et al.lll995t ) chose not to use an 
approximation for the binding energy of the envelope, but 
instead integrate the envelope mass from stellar structure 
calculations. As we will see in § 12.41 this choice does not 
necessarily lead to a more accurate value of a. So, here we 
determine as accurate an approximation to the binding en- 
ergy as possible. 

The binding energy of a star is: 



-Eta in 



Gm 



dm. 



(8) 



In the case of giants, the small and dense core, is surrounded 
by a vast envelope whose density rapidly falls with radius. 
Therefore we model the core as a point mass and the enve- 
lope as a shell of homogeneous density: 



r 47r(Ai?) 2 ' 

located at a distance XR from the core. The factor A < 



1 Different studies substitute the Roche Lobe radius for A{, in- 
stead of the actual separation between primary and secondary 
at the time of Roche lobe overflow. The initial orbital energy 
term is mostly negligible compared to the final orbital energy one 
though, so specific choices of the initial orbital separation are not 
fundamental. 



1 accounts for the fact that, in the shell approximation of 
the real envelope, the shell is located at the mass-averaged 
stellar radius. Using this assumption we calculate separately 
the envelope binding energy deriving from the attraction of 
the core, E cc , and that deriving from the attraction of the 
envelope onto itself, E ee : 

-Ebin = E cc + E cc . (9) 
The core-envelope binding energy can be easily calculated: 



E cc = — G 



M c dm _ M c M e 



XR 



XR ' 



(10) 



where dm is a parcel of envelope mass sitting at distance XR 
from the core, and the integration is over the entire envelope, 
approximated as a shell with radius XR. The self-gravity of 
the envelope results in: 



F - ° 
E cc - -- 



dm dm 



(11) 



where dm and dm are mass parcels in the envelope shell 
separated by a distance |r — r |. The factor | prevents dou- 
ble counting the interaction between parcels. Since we have 
made the assumption that all mass parcels are in a shell 
with radius XR, |r I = Irl = AT?, the denominator is: 



1 °° 



A^ 1 



(cos^O) i 



(12) 



where P n is the Legendre polynomial. For the final step, we 
only consider P = 1 and Pi = cos?/>, where ip is the angle 
between r and r and we neglect orders n > 1. Because of 
symmetry, the term including cos?/' vanishes. Substituting 
Eq. [12] into Eq. 111! we obtain the binding energy of the 
envelope in the shell approximation: 



Ebin ~ — G 



XR 



(13) 



This expression is different from all of the ones used 
before (the numerator of our expression is the lowest) al- 
th ough our express ion is actually quite similar to that used 
bv lWebbinkl || 19841 ). 



2.3 The thermal energy 

The virial theorem can be used to quantify another source of 
energy that may play a role in the CE interaction , namely, 
the thermal energy of the envelope (I WebbinklbOO&t ) . The fa- 
miliar identity 2K + U—0, where K is the total thermal, or 
kinetic energjjf] of the star and U is its potential (binding) 
energy, accurately represent the global properties of the en- 
tire star. If we only include the stellar envelope, rather than 
the entire star, the virial t heorem takes a s lightly different 
form (for the derivation see I Webbinkl |200l)): 



2K anv + ATiRt P c + E h 



0, 



(14) 



2 Throughout this paper we will keep referring to the energy due 
to the thermal motion of the particles as the thermal energy. This 
will avoid confusion with the energy stored i n the bonds of at oms 
and molecules, called the internal energy bv lWebbinkl fcOOSl '). 
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RGB 



AGB 



2t/ C nv (ergs) 
4itRIP c (ergs) 
2U eav + 47rfl3p c (ergs) 
-Ebin.onv (ergs) 
error without extra term 
error with extra term 



0.9983 X 10 48 
0.0532 x 10 48 
1.0516 x 10 48 
-1.0433 x 10 48 
4 % 
0.8 % 



2.7101 x 10 47 
0.1365 x 10 47 
2.8466 x 10 47 
-2.8348 x 10 
4 % 
0.4 % 



47 



Table 1. Error on the virial theorem calculations for the entire 
star using the 2 Mq model. R c defined by the left over mass 
criterion in § 12.41 



where K cnv is the thermal energy of the envelope only, and 
we now use -Ebm.cnv to represent the binding energy of the 
envelope only; R c and P c are the values of the radius and 
the pressure at the core-envelope boundary, respectively (de- 
fined in § [H| . 

We can show that the extra term (4-kRIPc) can be ne- 
glected, by integrating stellar structure models (whose de- 
tails will be given in § I2.4[) . We use a 2-Mq main sequence 
model evolved to the RGB and AGB phases. The total stel- 
lar binding energy is integrated using the denominator in 
Eq. 1161 A similar equation is used to determine the total 
thermal energy of the stellar models: 



U 



3 2 
, -P{r l )ATxr i An 



(15) 



where r» is the radius of the i th concentric shell, Ari are the 
shell's thickness, and P{r'i) is the pressure in the shell. R c is 
the core-envelope boundary radius chosen according to the 
criteria explained in § 12.41 

We see (Table [TJ that the first two terms in Eq. Q3] do 
equate to the third term to within less than 1% for both the 
RGB and AGB models. We also see that not including the 
second term, still results in an acceptable equality between 
twice the thermal energy and the gravitational binding en- 
ergy of the envelope. In this way we have a very simple way 
to account for the thermal energy of the stellar envelope in 
the CE energy budget: the thermal energy is simply one half 
of the binding energy. Since the two energy sources have op- 
posite signs, the effect of accounting for the thermal energy 
is to "fill" the stellar envelope potential well, or make the 
envelope lighter. 

Finally, one should remember that there is another pos- 
sible energy source: the dissociation a nd ionisation ener gy of 
the envelope. T his was discussed by lHan et all (|l995l ) and 
IWebbinkl (|200ST l. We will return to this topic in §H 



2.4 The stellar structure parameter A 

To determine a suitable value of A we use a 2-M model 
calculated with the cod e Modules for Exper iments in Stellar 
Astrophysics (MESA! 3 ! IPaxton et al.l l2010h . evolving from 
the main sequence to the tip of the AGB. This ne w stellar 
evolu tionary code compares well with EVOL (e.g. iHerwid 
120041 ). 

By equating the envelope binding energy obtained by 
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numerical integration of the stellar structure and that ob- 
tained by the approximation in Eq. 1131 we calculate A in this 
way: 

-{GM e /R)(M e /2 + M c ) 



A = 



-47rG Y^ n =R c m int{r l )p l r i Ar i 



(16) 



where pi is the density and radius of a shell of material of 
thickness Ar^ and m- llltl is the mass internal to that radius. 

The choice of a core radius based only on den- 
sity (e.g., where rfm /rfl ogr reaches a min imum, see also 
iTauris fc Dewill200ll and lBisscherouxlll99Sl ) locates a core- 
envelope boundary that is quite different from any boundary 
located on physical grounds. We therefore chose the core- 
envelope boundary to be at the radius where the nuclear 
burning reaches a maximum plus the thickness of a remain- 
ing envelope of mass 1CT 3 M for the 2.O-M AGB model 
described above and 10 -2 M for the corresponding RGB 
model. This criterion is motivated by the fact that AGB and 
RGB stars depart their respective giant branches by con- 
tracting, when the envelope mass decreases below a thresh- 
old value (or, f or RGB stars, when c ore helium is ignited; 
iBloeckerl [l995l ; ICastellani et alj 120061 ). In a CE situation, 
such contraction would detach the primary from its Roche 
lobe and dictate the end of the interaction. This criterion 
leads to very simila r core-envelope bound aries and values of 
A to the criteria of iDewi fc Taurisl (2000). For more discus- 
sion see Appendix lAl 

With this choice of core-envelope boundary, we carried 
out stellar evolution calculations for main sequence masses 
of 1, 1.5, 2, 2.5, 3, 4 and 5 Mq. For each of these models, we 
determined A for a few stellar structure models clustered on 
the RGB and AGB. The results are summarised in Table 2. 
For each main sequence mass, we list the numbers of stellar 
structures used in determining the average value of A, the 
smallest and largest stellar radius for the models in each 
cluster and finally the average value of lambda as well as 
the 1-cr spread of values. The minimum and maximum radii 
on the RGB and AGB for each model cluster were chosen as 
follows: the minimum radii were those of the first model that 
had developed a clear shell burning structure (hydrogen for 
the RGB and helium for the AGB). The maximum values 
were the maximum RGB and AGB radius values for the 
respective models. The value of A does not vary significantly 
within each cluster, despite the relatively large radius range, 
because the mass is centrally concentrated. 

We have fitted the average A values as a function of 
main sequence mass (Fig. [TJ and obtained the resulting fit 
parameters: 

/ M \ ( 2 - 11±al2 > 
Argb = (0.547±0.068)x ( ) 



A M © J 



xexp(-M MS /M )(17) 



and 



Aagb = (0.237 ± 0.021) + (0.032 ± 0.006) x 



Mms 
M„ • 



(18) 



which are valid for 1.0 M ^ Mms ^ 5.0 M . By calculating 
the value of A for each of our systems using the expression: 

A -1 ~ 3.000 - 3.816M e + 1.041Af 2 + 0.067M e 3 + 0.136M e 4 



of Wcbbink ( 20081 ) . we obtain much larger values. 

Finally, we comment on the fact that in our static cal- 
culation of the binding energy we assume that the enve- 
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# of models 


-^min 


-^m ax 


A 


(Mq) 


in cluster 


(Re) 


(R©) 




RGB 


1 


6 


20 


120 


0.175 ±0.055 


1.5 


6 


20 


120 


0.275 ±0.049 


2 


4 


20 


80 


0.352 ±0.030 


2.5 


5 


10 


30 


0.300 ± 0.028 


3 


4 


10 


40 


0.218 ±0.063 


4 


5 


10 


50 


0.146 ±0.032 


5 


4 


20 


80 


0.121 ±0.014 


AGB 


1 


4 


150 


300 


0.203 ± 0.073 


1.5 


6 


150 


400 


0.291 ±0.047 


2 


7 


50 


350 


0.302 ±0.052 


2.5 


6 


50 


300 


0.344 ± 0.039 


3 


7 


50 


350 


0.330 ± 0.047 


4 


9 


100 


500 


0.369 ± 0.047 


5 


8 


150 


500 


0.393 ± 0.027 



Table 2. Average values of A calculated for clusters of models 
equally distributed in radius between R m i n and i? max , on the 
RGB and AGB 

0.6 | , , , , , 1 

0.5 - 



0.4 - 




0.2 



0.1 

I 1 1 1 1 1 

1 2 3 4 5 6 

Figure 1. The average value of A (symbols with error bars) cal- 
culated for clusters of models on the RGB (top panel) and AGB 
(bottom panel) as a function of main sequence mass, along with 
the best fits (solid lines) 

lope that needs to be ejected is the same as the mass above 
the core-envelope boundary just discussed. This assumption 
does not take into account that in an evolving, mass-losing 
and expanding giant, the actual mass left on the star af- 
ter termination of the AGB may not be ex actly that which 
we ac counted for in our static calculation. iDelove fe Taaml 
|20ld ) calculate the adiabatic mass-loss response of the pri- 
mary. They explain that massive stars - they show the 
case of a 10 Mq - have a shallow entropy gradient in the 
core/envelope region, and this can lead to a 30% differ- 
ence between the donor's post-CE remnant mass and M c . 
However, for the lower mass stars we are interested in, the 
gradient is much steeper and two masses are similar. The 



other concern is that the binding energy as we calculate it, 
does not take into account the work done by the envelope 
throughout the adiabatic mass-loss. However, according to 
iGe et all (|20ld ). including this effect does not make a large 
difference on the binding energy (see their Section 4 and 
their equations 62 and 83). 



3 THE DETERMINATION OF a USING 
SIMULATIONS AND OBSERVATIONS 

Here we calculate a for observed post-CE binary systems as 
well as for a set of CE simulations. To do so, we use our ap- 
proximation (Eq. [5]) together with the value of A determined 
from the fits in § |2TH (Eqs. [171 and [TBI) . 

3.1 The pre-CE giant reconstruction technique 

We assume that observed post-CE systems with M c ^ 
0.47 M© have suffered a CE interaction on the AGB and 
those with M c < 0.47 Mq have suffered an interaction on the 
RGB. The latter assumption is not strictly correct, because 
more massive primaries can have a core mass >0.47 Mq 
during the RGB ascent (see for instance the 5-Mq model in 
Fig. [5} . However, our assumption is approximately correct 
because lower mass stars (< 2.5 Mq), which are statistically 
more common, have cores in the 0.47-0.62 Mq range only 
during their AGB ascent. In § 13.4.21 we discuss this topic 
further. 

Of the post-AGB primaries (M c ^ 0.47 Mq), some are 
surrounded by a PN, guaranteeing that the system only 
emerged from the CE interaction recently (or the PN, which 
has a lifetime of <100 000 years, would not be visible). 
When a PN is lacking, there is a chance that the time since 
the CE interaction might be sufficiently long to have al- 
lowed the evolution of the system's orbital period. In or- 
de r to avoid such systems we use the binary post-CE lis t 
of ISchreiber fc Gansickel l|2003l ) and lZorotovic et all (|201CJ ). 
who determined which post-CE systems have the same pe- 
riod today as they had when they emerged from the CE. 

Once the RGB or AGB origin of the primary is deter- 
mined, the radius and mass of the giant at the time of the 
CE interaction can be estimated. For post-RGB primaries, 
the radius has trad itionally been ob t ained by the core mass- 
radius relation of llben fc Tutukovl (|l985l . R ~ 10 3 5 M C 4 ). 
RGB masses at the time of the CE interaction cannot be 
determined fr om the core mass alon e, so a range has to be 
adopted (e.g., [Nclemans et al. 200^) 

To determine the pre-CE radius and mass for pri- 
maries that went thro ugh the interacti o n on the AGB, 
several studies, (e.g., lAfsar fc Ibanogiul l200ct l used the 
core mass-luminosity rel ation for post-AGB stars (e.g., 
IVassiliadis fc Wood! fl994l . L/L = 56694(M C /M - 0.5)) 
to determine the luminosity of the primary at the time of 
the CE interaction. They then used the initial-to-final mass 
relation to determine the main sequence progenitor mass of 
today's primary, to then determine the mass of the primary 
at th e time of the CE interaction (e.g.. IVassiliadis fc Woodl 
Il994h . Once the primary's mass and luminosity at the time 
of the CE interaction have been deter mined, the rad i us can 
be found using the fitting relation of iHurlev et al. I (|2000l . 
R = 1.125M-° 33 (L ' 4 + 0.383L 76 )). This method cannot 
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1 



0.9 




0.4 1 ' 1 1 1 1 1 

1 2 3 4 5 6 

M MS /M Sun 

Figure 2. A comparison of the MESA code initial-to- final 
mass relation (sol id line) and that determined empirically by 
IWeidemarml ||2000| . dashed line). 

be applied to primary stars with core mass < 0.55 Mg, 
be cause below that va l ue the core mass-luminosity relation 
of IVassiliadis fc Wood! (|l994l ) becomes extremely imprecise 
and for core masses ^ 0.5 M , the stellar luminosity be- 
comes negative. Finally, using the initial-to-final mass re- 
lation in this way is incorrect, because the CE interaction 
interrupted the evolution of the primary and the growth of 
its core, resulting in a less massive core (i.e., smaller final 
mass) than if the star had been single. 

By studying the relations used in the past one realises 
that reconstructing the giant's mass and radius using only 
the mass of the post-CE primary (assumed to have been the 
mass of the core of the giant), is fundamentally a statisti- 
cal process. Each relation appears to state that the error 
on a given quantity, e.g., the giant radius at the time of 
the CE interaction, depends solely on the error on the mea- 
sured core mass. However, this error propagation does not 
include the error in the fit, i.e., the error in the relations 
themselves. Nor is the error in the assumptions used to es- 
tablish the relations quantified. We therefore diverge from 
previous methods, and recalculate the fitting relations from 
stellar evolution models. This allows us to better quantify 
the true uncertainty on the entire reconstruction process. In 
our reconstruction technique, we make use of two sets of stel- 
lar evolution calculations. The first set is from the detailed 
calculations of MESA for Z = 0.01 (§|2~3|) The second set is 
from the computations of iBertelli et al.l l|2008l ). who calcu- 
lated an extensive grid of models including their metallicity 
dependence. 

3.1.1 The primary mass at the time of the CE interaction 

To determine the mass of the primary at the time of the 
CE interaction for the post-AGB group, we start by deter- 
mining the mass of the main sequence progenitors. To do so 
we use the initial-to-final mass relation (IFMR) determined 
with the MESA code. The final mass is assumed to be the 
mass of the core at the first thermal pulse. Although this 
may appear as an under estimate of the final mass, it is ac- 
curate for core mass values < 0.52 Mq, because low mass 
stars with low mass cores have low mass envelopes which 
get depleted rapidly by mass loss and do not give the core 
time to grow. Our assumption also results in good accu- 
racy for masses > 0.58 M®, because high mass stars suffer 



substantial dredge-up so that the core does not grow appre- 
ciably during thermal pulses. For values between 0.52 and 
0.58 Mq, the core could grow slightly above the value that 
we are using, resulting in a slight overestimate of the initial 
mass. However this uncertainty is very small compared to 
the very large uncertainty on the core mass, as we explain 
below. 

Older predictions fo r the IFMR (e.g., iBloeckerl Il995l ; 
IVassiliadis fc Wood|[l99l ) compare well to the MESA IFMR 
for lower initial masses (<3 Mq), while for the larger ini- 
tial masses the more rece nt MESA IFMR y ields higher WD 
masses (also observed by lHan et al.l |l994)). The reason is 
that the new MESA models include convective core over- 
shooting during the main sequence, which increases the core 
mass predictions compared to the older models. However, 
at low initial mass, models now show efficient third dredge- 
up, which compensates for the core mass increase from core 
overshooting. For more massive AGB stars the third dredge- 
up does not play a significant role in the core mass predic- 
tion, and the core mass difference from core overshooting 
assumptions can be observed. The MESA IF MR compares 
well w ith the observationally derived IFMR of I Weidemannl 
( 2000, Fig. 0, although the observational IFMR has a large 
scatter. 

For each system's primary we obtained the main se- 
quence mass of its progenitor by direct interpolation. The 
metallicity dep endence of the IFMR is small, as can be seen 
by a fit to the IBertelli et al.l (|2008h models (M W d/M q = 
0.451 + 0.091Mms/Mq - 0.044 log Z/Z e ) and we ignore it. 
Before we can use the IFMR, we need to account for the fact 
that the IFMR relates the final WD mass, Mwd, to its main 
sequence mass, Mms, only if the star evolved to the natural 
termination of its AGB life. The CE interaction interrupted 
the regular AGB evolution and core growth, so that the 
post-CE primary mass, M c , is smaller than it would have 
been had the primary been si ngle (Mwd). This m ass dis- 
crepancy is determined for each IBertelli et aD l|2008h model, 
assuming that a CE will terminate the AGB evolution at 
a random value of R, between the m aximum RGB radius 
and the maximum AGB radius. We use IBertelli et all (<2008l ) 
models in the mass range 0.8-2.5 Mq, and metallicity range 
Z = 10" 4 - 0.07. On average, models depart the AGB with 
a mass that is 0.028±0.024 M lower than if the AGB evo- 
lution had progressed till its natural termination. The error 
in this value is dominated by the lack of knowledge of the 
precise stage at which the system departs the AGB. 

We can now calculate the WD mass of each of our post- 
AGB, post-CE primaries and use the IFMR to determine 
the main sequence mass of its progenitor. The error on the 
WD masses is determined from the error on the measured 
post-CE primary masses and the error on the core growth 
estimation, added in quadrature. In Fig. [3] (left) we show the 
MESA IFMR with our interpolated data points and their 
error bars. To determine the error on the main sequence 
masses thus obtained, we projected the errors calculated on 
Mwd on the IFMR itself, as demonstrated in Fig. [3] (right). 
For three data points at the higher (V471 Tau) and lower 
(SDSS J0110+1326 and RR Cae) extremes of the IFMR, the 
errors were assumed to be the same as the relative error of 
HS 0705+6700 and SDSS J1548+4057, for the low and high 
mass, respectively, because these two systems have the same 
measurement errors and similar core mass values. Using this 
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Figure 3. Left: the initial-to-final mass relation (IFMR) determined from the MESA code calculations, plotted together with the data 
and their error bars (Table[3]l. The dotted lines delimit the estimated error on the IFMR itself. Right: the IFMR plotted with the data 
point for BE UMa to exemplify our determination of the error on the derived main sequence masses. 



error estimation method, we are not taking into account 
the error on the IFMR itself (plotted as a dotted curve in 
Fig. [3l from approximate estimates). If we had, the errors 
on the determined main sequence masses would have been 
only marginally larger. 

The last step in the main sequence mass determination 
is carried out by using conditional probability considera- 
tions. For each of the three populations of post-CE binaries 
(the central stars of PN, the post-AGB stars with no PN 
and the post-RGB stars) we know the main sequence mass 
distribution of the progenitors from population considera- 
tions and observations (see below). To exploit this extra 
knowledge, we convert the progenitor main sequence mass 
of each system's primary, determined above from the WD 
mass alone (dotted vertical lines in Fig. [4]), into a Gaussian 
probability distribution function (PDF), by using the error 
estimates (dotted curves in Fig. 2}. We then multiply each 
PDF with the main sequence mass distribution for the par- 
ent population (also called the prior in Bayesian statistics; 
dashed curves in Fig. This procedure results in a PDF for 
the "conditional" values of the main sequence masses of each 
of our post-CE primaries (solid curves in Fig. 3]), from which 
we can determine a mean value (solid vertical lines in Fig. 2]) 
and a new error estimate. These are the values we adopt for 
the primary progenitors' main sequence masses. We note 
that two identical primary masses (M c ) that have differ- 
ent error estimates will result in two different "conditional" 
main sequence progenitor masses, because of the Bayesian 
statistical treatment (compare HW Vir, HS0705+6700 and 
MS Peg in Table[3]). In Fig. [4] we show two specific cases, one 
for a central stars of PN, and one for a post-AGB system 
with no PN. 

In Fig. [4] (dashed lines) the prior distributions are the 
main sequence mass distribution of the parent populations. 
For the central stars of PN, the progenitor main sequence 
mass distributio n is known from populati on synthesis as well 
as observations (|Moe fc De Marco] I2006I . dashed line in the 
left panel of Fig.QJ. For the post-AGB stars with no PN, the 
main sequence mass distribution of the proge nitors can be 
deter mined using the WD mass distribution (| Kepler et al.l 
l2007h and using the IFMR to translate the WD mass dis- 
tribution into a main sequence mass distribution. The WD 



mass distribution of iKepl er et al. is fitted by three 

Gaussian curves with different means and standard devia- 
tions. The Gaussian curve with the lowest mean mass cor- 
responds to the helium WDs and we ignore it. To represent 
the WD mass distribution corresponding to the post-AGB 
stars (with no PN), we only use the remaining two Gaussian 
curves combined (with mean masses of 0.578 and 0.678 Mq, 
respectively, standard deviations 0.047 and 0.148 Mq, re- 
spectively and relative strength 2.9: lj3- 

For the post-RGB group we cannot use the IFMR, be- 
cause the primaries have not been through the AGB evo- 
lution. For all post-RGB primar ies, we theref ore use the 



initial mass function (oc M~ 2 35 : iKroupal 1200 ill truncated 



between 0.78 and 2.3 M (|Nelemans et al'l feoooi This dis- 
tribution results in a mean value accompanied by an error 
(1.19±0.40 M©). 

Once we have calculated the primaries' progenitors 
main sequence masses, we can derive the mass of the primary 
at the time of the CE interaction (Mi), by determining how 
much mass the star lost between the main sequence and the 
time of the CE interaction on the RGB or AGB . The stellar 
evolutionary calculations of lBertelli et~aH |20Qg| ) for a range 
of masses and metallicities are used once again to average 
the RGB and AGB stellar masses using uniform weighting 
with respect to radius (for the AGB, we used only the values 
for which the radius was larger than the maximum radius at- 
tained on the RGB). The values thus obtained are Mi/Mms 
= 0.98±0.02 for post-RGB stars and 0.87±0.07 for post- 
AGB stars, where the errors are derived from the scatter on 
the fit. We then considered that the presence of a companion 
would stimulate ma ss-loss prior to the CE interaction (e.g., 
iBear fc Soken feoiO), so we lowered these numbers slightly 



4 We note that this prior is not entirely correct, because CE inter- 
actions on the AGB happen slightly more frequently to relatively 
more massive stars. This is because lower mass stars grow to rel- 
atively larger radii while on the RGB, increasing the chance that 
an interaction will happen on the RGB rather than the AGB. As 
a result the priors for the post-AGB, post-CE binary populations 
should be slightly weighted towards more massive stars. To calcu- 
late the correct priors would entail a complete population study, 
which we have not carried out here. 
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Figure 4. The main sequence mass (solid vertical lines) determined for two of our post-CE primaries: A 63, a central star of PN (left) 
and V471 Tau, a post-AGB star (right). These masses are the means derived from statistical distributions (solid curves), that are the 
result of a convolution of the Gaussian mass obtained from the IFMR (dotted curve with the mean plotted as dotted vertical lines) and 
the Baysian prior distributions (dashed lines) . These latter distributions are the mass distributions of the progenitors of the two types 
of objects, the central stars of PN (left) and the post-AGB stars with no PN (right). 



and raised the uncertainty: Mi/M M s = 0.90±0.10 for RGB 
systems and 0.75 ± 0.15 for AGB ones. Once again, the de- 
pendence on the mass and metallicity is not what dominates 
the error, but rather the lack of knowledge of what stage dur- 
ing the RGB or AGB evolution the systems enter the CE 
phase. 



3.1.2 The primary radius at the time of the CE 
interaction 

We next determined the radii of the giants at the time of the 
CE interaction. From the MESA calculations we know the 
radius evolution as a function of the mass of the hydrogen- 
exhausted core (Fig. [5]). We first selected the track corre- 
sponding to the main sequence progenitor's mass of each of 
our post-CE primary stars (§ 13 . X . 1 p : we then read from the 
plot the radius corresponding to today's core mass. This is 
the radius of the giant at the time of the CE interaction. 
The error on this radius has to be determined taking into 
account the errors for both the measured core mass and the 
derived main sequence mass. Considering each error bar in 
turn results in new radius values, from which we can esti- 
mate the final error on the radius. There is, however, one 
additional complication. The radius evolution of each model 
is complex and non-monotonic. This means that we need to 
use logical arguments to consider each radius value obtained 
in this way: a star can only be caught in a CE interaction 
while its radius is growing and if it is larger than at any 
time in the past. Once all these conditions are accounted 
for, it is found that the range of possible radii correspond- 
ing to each of our systems is quite large and that it is no less 
acc urate to use fi t s of th e core mass vs. radius values from 
the iBertelli et al.l ((2008) calculations, for the usual ranges 
of masses and metallicities: 

R = 44Oi? Q (M M s/M Q )- o ' 47 (A/ c /(O.6A/ )) 5 ' 1 (Z/Z ) - 15 , (19) 

where the rms error on the fit is ±50 R . Or, evaluating R 
in terms of Mi instead of Mms- 



R = 440i? Q (M 1 /M Q )- a54 (M c /(0.6M Q )) 5 '°(Z/^ Q ) ai4 , (20) 

with the same fit scatter. We note that this relation applies 
equally to the post-RGB and post-AGB stars, although the 
fit was technically carried out for the AGB phase only. These 
fits are over-plotted on the detailed MESA calculation in 
Fig. [S] As can be seen from Fig. [5] the radius values implied 
by our approximate fits and the range of radii one would 
derive using the detailed MESA tracks become quite differ- 
ent for larger main sequence masses. The most massive of 
our systems, V 471 Tau, with a core mass of 0.84 Mq, has a 
conditional main sequence mass of 3.9 Mg. Using the MESA 
tracks one derives a radius at the time of the CE interaction 
of 435T100 R0. The value derived instead from the fits is 
570±62 (see Table [3]). These two values are quite different, 
although they are consistent within the uncertainties. Pre- 
emptying our results from § 13.41 the MESA radius value re- 
sults in a=0.28±0.16, while using the radius fit from Eg. 1201 
a=0.21±0.15. We therefore see that even for the worse case 
scenario of the system with the largest mass, the advantage 
of using the detailed MESA code for this purpose is limited. 

Finally, as we did for the case of the main sequence 
mass, we combine these radius values and their formal er- 
rors with the prior knowledge that for our primaries' progen- 
itor populations, the radii of the RGB primaries should be 
in the range 10-300 Rq, while for t he AGB primaries t hey 
should be in the range 50-650 R ^Bertem et al.ll2008T ), in 
the usual mass and metallicity ranges. The maximum RGB 
radius of 300 Rq , and maximum AGB radius of 650 Rq , are 
the maximum radii achieved by any model in the set con- 
sidered. The minimum radius value adopted for the AGB 
(50 Rq) was selected because it is the smallest of all maxi- 
mum RGB radii, achieved by any of the considered models. 
The minimum RGB radius (10 Rq) is the approximate ra- 
diu s that separat e s Her tzsprung gap from RGB models for 
the lBertelli et afl (|2008f ) tracks. Convolution with the prior 
results in a new radius value and a new error, according to 
Bayesian statistics. 
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3.2 Observed systems used in the 
determination of a 



3.3 The simulations used in the 
determination of a 



In Table [3] we list the values of M c , M2, P and Af deter- 
mined from observations as well as the values of Mi and R 
deduced from stellar evolution in § 13.1) References for the 
observed values can be obtained from iDe Marcd HH), for 
central stars of PN. Of al l the ~ 40 post-CE binary central 
stars listed by IDe Marcd |2009l ) only 6 central stars with 
non-degenerate secondaries have sufficient information to be 
useful in the present study. For the observed post-CE binary 
in the centre of PN NGC 6337 we list two sets of parameters 
in Tables [3] and H When we fit the value of a (§ l3A2j) . we 
use only one set (that with M2 = 0.35 Mq). The fit with the 
alternative set, results in very similar fit parameters. Finally, 
for the binary syst em in the PN H FG 1 there also are two 
sets of parameters |De Marcdl20091 ). We use the set with the 
most massive companion, M 2 = 1.09 Mq, becau se a recent 
analysis of X-ray radiation l|Montez et al.ll2010h finds this 
system's companion to be coronally active and the only way 
to reproduce the X-ray luminosity is to have a companion 
at least as massive as a ~ 1 Mq main sequence star. 

Other post-CE bina r ies ar e primarily from the list 
of ISchreiber fc Gansickd l|2003l . where we selected only 
those systems that have not suffered substantial reduc- 
tion of the orbital period after emerging from the CE). 
Whereve r parameters have been updated by new measure- 
ments of IZorotovic et al.l l|2010l ) , we have adopted the re- 
visions. From the extensi ve list of Sloan Digital Sky Sur- 
vey post-CE binaries in IZorotovic et al.l ^20 10), we only 
use the systems SDSS J1548+4057, SDSS J0110+1326 and 



SDSS J 1435+3733 (for more disc ussion, 



\B- F inally, 



we have also considered AA D or (|Miiller et al.l I2OI0I ) and 
HD 149382 (|Geier et al.ll2009l ). because of their low mass 
secondaries. Whenever errors were not given in the cited 
references, we have assumed an error of 20% on both pri- 
mary and secondary masses. In Table[4]we list the calculated 
values of q (—M2/M1) and a, determined using Eq.[6] along 
with their uncertainties, both in linear and l ogarithmic form . 
We also list values of a determined using the lWebb ink (2008) 
relation. 



Hydrodynamic simulations, carried out both with smooth 
particle hydro dynamics (e.g. Terman Sz TaarrJl 19961 ) or grid 
methods fe.g.. lSandouist et al.lll998h . are atool to determine 
the value of a ab initio. These simulations start by mapping 
a giant star into the simulation domain. The quantities de- 
scribing the star (i.e., density, temperature, internal energy, 
etc.) are calculated using a one-dimensional stellar evolu- 
tion model. The companion is represented by a point mass, 
as is the compact core of the giant. The physics included 
in the simulation is typically only the gravitational attrac- 
tion of all the masses involved. The interaction timescale 
is short enough that radiation transport can be neglected 
(e.g., see discussion in Sandquist et al. 1998). The in-spiral 
of the companion is typically followed until the resolution of 
the grid becomes insufficient, which may or may not be at 
the termination of the CE interaction (i.e., when sufficient 
mass has become unbound from the system to dictate the 
collapse of the primary and its detachment from its Roche 
lobe). At the end of a simulation, a can be determined by 
using Eq. [6] 

We determ i ned a from 5 simulati ons, four from 
ISandquist et all (|l998l ) and one from Marco et aJj 
J2003t). All simulations were carried out with the 
Burkert fc Bodenheimeil d 19931 ) code, as modified by 



Sandquist et al.l (| 19981 ). Relevant parameters are taken from 



the list ed publications and are reproduced in Tables[3]where, 
for the lSandquist et al.l jl998) simulations, we use the same 
labels as in their Table 1. Their simulations 1 and 2 were 
identical, except for the giant envelope's rotation being ei- 
ther zero or synchronised with the companion's orbital mo- 
tion. For these two simulations we take the average of their 
post-CE binary periods. All simulations considered were 
stopped after orbital decay time-scales lengthened consider- 
ably, when only a very small amount of envelope mass was 
left inside the orbit of the companion (~ 10 -4 - 10 -3 Mq). 
However, in all cases a considerable amount of envelope gas 
was still bound to the system albeit at some distance from 
the core (usually outside the original giant radius and with 
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LM Com 
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0.45± 


0.05 


0.28± 0.05 


0.26 


1.5± 0.1 


1.2± 0.4 


1.1± 


0.4 


116zfc 61 


RR Cae 


pRGB 


0.44± 


0.02 


0.18± 0.01 


0.30 


1.6± 0.1 


1.2± 0.4 


l.lzfc 


0.4 


89zfc 29 


HD 149382 


pRGB 


0.47± 


0.12 


0.015± 0.007 


2.39 


5.9± 1.5 


1.2± 0.4 


l.lzfc 


0.4 


136zfc 91 


HR Cam 


pRGB 


0.41± 


0.01 


0.10± 0.01 


0.10 


0.72± 0.02 


1.2± 0.4 


l.lzfc 


0.4 


59zfc 13 


CC Cet 


pRGB 


0.39± 


0.10 


0.18± 0.05 


0.28 


1.5± 0.3 


1.2± 0.4 


l.lzfc 


0.4 


96zfc 82 


WD 0137-349 


pRGB 


0.39± 


0.04 


0.052± 0.005 


0.08 


0.60± 0.05 


1.2± 0.4 


l.lzfc 


0.1 


54zfc 26 


DeMa 2 


Sim 


0.60 


0.10 


30 


36 


1.20 


1.0 


645 


Sandl-2 2 


Sim 


0.70 


0.40 


1.08 


4.56 


3.18 


2.S 


190 


Sand3 2 


Sim 


1.00 


0.40 


0.90 


4.39 


5.29 


4.f 




190 


Sand4 2 


Sim 


1.00 


0.60 


0.96 


4.79 


5.29 


l.f 




190 


Sand5 2 


Sim 


0.94 


0.60 


2.48 


8.90 


5.31 


l.( 




353 



1 CSPN: the primary is a central star of PN and went through the AGB evolution; pAGB: the primary 

does not have a PN, but we consider it as a pot-AGB star; pRGB: the primary is considered as having suffered a 

CE interaction on the RGB and never having ascended the AGB; Sim: simulations. 

2 DeMa: Dc Marco ct al. (2003). Sand: Sandquist ct aL (1998), simulations. The number indicates the particular simulation, 
and follows the scheme of their Table 1 . See also text in § 13.31 



Table 3. The observationally-derived quantities for our post-CE systems (M c , Mi and P) listed alongside the quantities derived from 
orbital, stellar evolution and population considerations (j4/, Mi and R). Simulation inputs (M c , M2, Mi and R) and outputs (P and 
Af) are also listed. All these quantities are inputs to the a equation (Eq. [SJ> . Parameters in brackets are an alternative set not used in 
the fits (see § EH) 



small binding energy). Because of this left over bound enve- 
lope, we may not assume that the separation of the binary at 
the end of the simulation is the actual post-CE binary sepa- 
ration. Any envelope still bound to the system at the end of 
the CE will fall back onto the system form a circumbinary 
disk that might have some dynamical effect on the binary 
period. To understand this eventuality, will require further 
numerical simulations. 

The simulations of Yor ke et al. (1995) started the in- 
teraction with the companion well inside the AGB stel- 
lar envelope. This results in a negative value of a because 
M c /Af < (M c + M e )/A l . Although they justified this initial 
setup by noting that spin-up only becomes important within 



their initial radius, when the timescale for orbital decay ex- 
ceeds the orbital period, they also stated that a value for a 
cannot be trivially derived from their simulations, but that 
the values can be estimated between 0.3 and 0.6. Since it is 
not clear how these values were determined, we did not use 
their models. 

Finally, we did not use the simulations of 
iTerman fc Taaml (Il996l ). They considered a 5 M AGB 
giant in two different evolutionary stages. The less evolved 
AGB star, still at the very base of the AGB, and still 
undergoing core helium burning, was used to simulate a 
CE interaction with a 0.5 Mq companion and presumably 
resulted in a merger. The more evolved of the two stars, was 
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Name 


Type 1 


logg 


log a 


q 


a 


a(Web) 


BE UMa 


CSPN 


-0.85+ 0.12 


0.08+ 0.25 


0.14+ 0.04 


1.2+ 0.7 


1.1 


A 46 


CSPN 


-0.78± 0.13 


-0.45+ 0.34 


0.17+ 0.05 


0.35+ 0.31 


0.24 


A 63 


CSPN 


-0.78± 0.17 


-0.54+ 0.34 


0.16+ 0.06 


0.29+ 0.25 


0.23 


HFG 1 


CSPN 


0.01+ 0.16 


-1.06+ 0.35 


1.0+ 0.4 


0.087+ 0.080 


0.043 


DS 1 


CSPN 


-0.99± 0.10 


-0.29+ 0.23 


0.10+ 0.02 


0.51+ 0.28 


0.46 


NGC 6337 


CSPN 


-0.59+ 0.18 


-1.01+ 0.36 


0.26+ 0.11 


0.10+ 0.09 


0.068 






(-0.83± 0.19) 


(-0.83+0.34) 


(0.15+0.06) 


(0.15+0.14) 


(0.11) 


V471 Tau 


pAGB 


-0.50+ 0.15 


-0.68+ 0.30 


0.32+ 0.11 


0.21+ 0.15 


0.17 


UZ Sex 


pAGB 


-0.91+ 0.20 


-0.24+ 0.45 


0.12+ 0.06 


0.57+ 0.70 


0.37 


SDSS J1548+4057 


pAGB 


-1.13± 0.10 


-0.32+ 0.22 


0.074+ 0.018 


0.48+ 0.26 


0.31 


RE J1016-053 


pAGB 


-1.11+ 0.12 


0.13+ 0.25 


0.078+ 0.023 


1.3+ 0.8 


0.85 


Feige 24 


pAGB 


-0.56± 0.14 


0.08+ 0.29 


0.27+ 0.09 


1.2+ 0.9 


0.86 


IN CMa 


pAGB 


-0.52± 0.14 


-0.32+ 0.30 


0.30+ 0.10 


0.47+ 0.35 


0.35 


RE J2013+400 


pAGB 


-0.88± 0.16 


-0.16+ 0.31 


0.13+ 0.05 


0.69+ 0.53 


0.44 


NN Ser 


pAGB 


-0.94± 0.15 


-0.54+ 0.25 


0.11+ 0.04 


0.29+ 0.18 


0.17 


GK Vir 


pAGB 


-0.98± 0.13 


-0.22+ 0.30 


0.10+ 0.03 


0.60+ 0.45 


0.34 


Hz 9 


pAGB 


-0.65± 0.17 


-0.40+ 0.41 


0.22+ 0.09 


0.40+ 0.43 


0.28 


AA Dor 


pAGB 


-1.17± 0.20 


-0.16+ 0.44 


0.07+ 0.03 


0.69+ 0.81 


0.40 


BPM 6502 


pAGB 


-0.72+ 0.13 


-0.47+ 0.33 


0.19+ 0.06 


0.34+ 0.28 


0.21 


PG 1017-086 


pAGB 


-1.06+ 0.14 


-0.60+ 0.33 


0.087+ 0.027 


0.25+ 0.21 


0.14 


NY Vir 


pAGB 


-0.78+ 0.14 


-0.78+ 0.33 


0.17+ 0.06 


0.17+ 0.14 


0.10 


SDSS J1435+3733 


pAGB 


-0.59± 0.13 


-0.87+ 0.27 


0.26+ 0.08 


0.13+ 0.09 


0.084 


HW Vir 


pAGB 


-0.81± 0.16 


-0.73+ 0.40 


0.16+ 0.06 


0.19+ 0.20 


0.11 


HS 0705+6700 


pAGB 


-0.76± 0.15 


-0.82+ 0.34 


0.18+ 0.06 


0.15+ 0.13 


0.087 


MS Peg 


pAGB 


-0.48± 0.11 


-0.91+ 0.27 


0.33+ 0.08 


0.12+ 0.08 


0.074 


SDSS J0110+1326 


pAGB 


-0.31± 0.12 


-0.82+ 0.28 


0.49+ 0.13 


0.15+ 0.11 


0.10 


LM Com 


pRGB 


-0 58+ 15 


-0 43+ 34 


26+ 09 


37+ 32 


0.26 


RR Cae 


pRGB 


-0.77+ 0.14 


-0.09+ 0.29 


0.17+ 0.05 


0.81+ 0.57 


0.54 


HD 149382 


pRGB 


-1.85+ 0.21 


1.52+ 0.47 


0.014+ 0.007 


33+ 43 


16 


HR Cam 


pRGB 


-1.03+ 0.14 


0.04+ 0.27 


0.093+ 0.030 


1.1+ 0.7 


0.68 


CC Cet 


pRGB 


-0.77+ 0.17 


-0.09+ 0.45 


0.17+ 0.07 


0.82+ 1.02 


0.56 


WD 0137-349 


pRGB 


-1.31+ 0.14 


0.30+ 0.33 


0.049+ 0.016 


2.0+ 1.7 


1.2 


DeMa 2 


Sim 


-1.02 


0.53 


0.096 


3.4 


1.8 


Sandl-2 2 


Sim 


-0.84 


0.34 


0.14 


2.2 


1.6 


Sand3 2 


Sim 


-1.06 


0.96 


0.087 


9.0 


6.3 


Sand4 2 


Sim 


-0.88 


0.82 


0.13 


6.6 


4.8 


Sand5 2 


Sim 


-0.89 


0.86 


0.13 


7.2 


5.4 


1,2 : see comments to Tabic 3 



Table 4. The determined values of a, using our equation as well as that of IWebbinkl |l984, 2008). Parameters in brackets are from an 
alternative sets of inputs (see Table [3] and § 13.21 1 . 



used to simulate CE interactions with 0.5 and 1.0 Mq main 
sequence companions. In neither case was the simulation 
followed till the end of the interaction, and a substantial 
amount of envelope remained bound to the system at the 
end of the simulation. However, the orbital decay time-scale 
had considerably slowed down and envelope mass was 
still being unbound. The authors therefore concluded that 
the system would survive, albeit with a smaller orbital 
separation than that reached by the end of the simulation. 



3.4 Results 

Here, we finally calculate the values of a from observations 
and simulations (Table 3] where we report both linear and 
logarithmic values). The error bars on a are significant but 
realistic. 



3.4.1 a > 1. 

For 6 observed systems and all the simulations a exceeds 
unity. The physical interpretation of this is that the specific 
CE interaction under consideration has benefited from more 
energy than the orbit al energy supplied by the comp anion. 
Several authors (e.g.. lHan et al.lll995l ; [Webbinkll2008l ) dis- 
cuss thermal energy as well as the dissociation and ionisation 
energies as additional sources of energy in the a equation. In 
§ 12.31 we discussed the thermal energy and showed that we 
can use the virial theorem to gauge its value, which is exactly 
half of the binding energy. If we include the thermal energy 
in this way, all values of a will be below unity, except for the 
system with a substellar companion (HD 149382, discussed 
further in §2| and the simulations. 



3.4-2 A possible \ogq vs. log a anti- correlation 

In Fig. [6] we explore possible correlations between a and 
other relevant parameters. The simulations' a values are 
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Figure 6. A plot of log a as a function of log Mi (top), logM2 
(middle) and logarithm of the period (bottom). The symbols are 
the same as in Fig.[7J Error bars were omitted for clarity but were 
used in the linear fits (dashed lines). Post- RGB systems were not 
fitted in the upper panel. Simulations (crosses) were never fitted. 

plotted (cross symbols) but not fitted, because it is impossi- 
ble to determine their error bars in a way that is consistent 
with the observations. We fit log Mi vs. logcv for post-AGB 
and CSPN systems (not for post-RGB primaries, because for 
those stars we had to adopt a single mass value). We find 
a very marginal correlation (correlation coefficient, r=0.56). 
On the logM2 vs. logo plane we see a slightly better anti- 
correlation with r = —0.79 (post-RGB systems are fitted 
here). Fits to the data and their errors result in reduced 
X 2 ~0-77 and 1.79 for the Mi and M2 cases, respectively. 
The correlation between logo and the period is extremely 
weak (r = 0.33). 



A fit to logo as a function of log(q = M%jM\) is plot- 
ted in Fig. [7] In Table [5] we list the parameters of the linear 
fit with their errors, along with the number of data points, 
reduced \ 2 value, probability to exceecjf] and the correlation 
coefficient r; we then list constant a fits, along with their 
reduced x 2 an d probability to exceed. First, we fit all of 
our data together, except for the simulations, as before. The 
data is reasonably anti-correlated. We then exclude the out- 
lier (HD 149382 [q = 0.014, a= 33]); the anti-correlation is 
somewhat degraded, although these two fits are very equiv- 
alent from a statistical point of view. Both have a reduced 
X 2 below unity and a probability to exceed of 0.8, meaning 
that the errors may have been overestimated or too many 
degrees of freedom have been considered in the fit. How- 
ever, a constant a fit never represents the data better with 
a probability to exceed of zero, which means that random 
deviations would always explain the data better than the 
constant a fit. 

To investigate further our data, we eliminated groups 
of systems to determine whether this would improve the 
statistics of the anti-correlation among the remaining data. 
By fitting only post-RGB and central stars of PN systems 
can one improve the anti-correlation (Fig. [7J — right panel). 
On the other hand the probability to exceed reduces further 
and since the constant a fit is not better, we once again have 
to conclude that the errors have been overestimated for this 
fit. The improved anti-correlation might derive from the fact 
that in our reconstruction technique we have assumed that 
any primary with mass larger than 0.47 Mq, went through 
the AGB evolution. That assumption results in a specific 
primary mass at the time of the CE interaction (§ I3.1.1|) . It 
is however possible, that some of these primaries are actually 
post-RGB objects descending from more massive stars, as 
we described in § 13.11 More massive stars develop a helium 
core more massive than 0.47 Mq before ascending the RGB. 
Once they expand, they may suffer a CE interaction on the 
RGB that truncates the RGB evolution by ejecting the en- 
velope. These massive, post-RGB, post-CE primaries would 
later ignite core helium, thus developing a CO core, but with 
an envelope mass low enough to make them appear as hot 
sub-dwarfs with a close companion. It is therefore possible 
that these objects' masses and radii at the time of the CE 
interaction might have been significantly miscalculated, re- 
sulting in the observed large scatter on the logg vs. log a 
plot in Fig. [7] (left panel). On the other hand, the central 
stars of PN are (almost) guaranteed to have been on the 
AGB, and the systems with M c < 0.47 M Q are post-RGB 
stars, making these two groups less prone to erro r. 

The values of a derived by ISandauist et "all (| 1998T ) are 
approximately a factor of fo ur lower than those we derive . 
They used the formalism of ITutukov k, Yungelsonl (|l979l ). 
except that they revised the expression for the orbital en- 
ergy to include the orbital energy value at the t ime o f 
the CE interac t ion, u nlike ITutukov fc Yungelsonl (|l979T) . 
ISandauist et al.l (| 19981 ) also derived the binding energy as 

5 The probability to exceed is the probability that a set of ran- 
domly sampled data from a normal distribution produces a re- 
duced x 2 that exceeds the one determined. This value should be 
as close as possible to 0.5, corresponding to a reduced \ 2 °f a P" 
proximately unity, wh ere the exact value depends o n the number 
of degrees of freedom (Bevington & Robinson 2003). 
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Linear fit Constant fit 





71° 


r a 




Iojj 


; o= 


x> 6 


pte° 


a = 


x 2 /» c 


pte° 


All data 


31 


-0.78 


(- 


-1.4 ±0.2) - 


(1.2 ±0.2) logg 


0.78 


0.79 


0.43±0.08 


1.9 





As above, no HD 149382 


30 


-0.64 


(- 


-1.2 ±0.2) - 


(1.0 ±0.3) logg 


0.75 


0.83 


0.41±0.06 


1.4 


0.06 


Only CSPN and pRGB 


12 


-0.92 


(- 


-1.4 ±0.3) - 


(1.4 ±0.3) logg 


0.54 


0.86 


0.63±0.12 


2.8 





As above, no HD 149382 


11 


-0.83 


(- 


-1.2 ±0.3) - 


(1.2 ±0.4) logg 


0.53 


0.85 


0.54±0.11 


1.7 






a n: size of the dataset. r: correlation coefficient, pte: probability to exceed. 
'Degrees of freedom, v = n — 2. 
c Degrees of freedom, v = n — 1. 



Table 5. Statistical properties of the fits to logg vs. log a. 



the difference between the binding energy at the beginning 
and at the end of the CE interaction, where they integrated 
the stellar structure for the initial value and used the simula- 
tion output to determine the remaining binding energy from 
the left over bound envelope in the system. Although we sus- 
pect that the leftover envelope binding energy is low, they 
do not report the actual value. It is not possible at this point 
to account exactly for this discrepanc y. However a large part 
of it i s cert ainly due to the f act th at pTutukov fc Yungelsonl 
i| 19791) . like lYungelson et all (Il993p (Eqs.|2]and|5) use a far 
larger radius of the primary at the time of the CE interac- 
tion, reducing the value of a (see § I2.1|l . 



4 THE STELLAR RESPONSE AND THE 
THERMAL ENERGY 

In § !2.3l we discussed using the virial theorem to quantify the 
entire energy budget of the envelope. If we accounted for the 
thermal energy of the envelope in this way, we would mul- 
tiply the envelope binding energy by a factor of 0.5, which 
would lower all a values by a factor of two, and bring all but 
one of the observed systems below unity. 

It is possible that the anti-correlation in the data, where 
some of the systems with low q value appear to need the ther- 
mal energy, while others do not, may hold some clues. Since 
the virial theorem is based on the principle of hydrostatic 
equilibrium, and since adjustments to hydrostatic equilib- 
rium take place on the stellar dynamical timescale, we may 
hypothesise that the timescales of envelope penetration and 
in-spiral vary as a function of q. 

The equation of hydrostatic equilibrium dictates that, 
during their giant phases, stars expand upon losing mass 
because the pressure gradient overcomes gravity. When the 
envelope mass is reduced below a certain threshold by mass- 
loss, the gravity term dominates and the star shrinks rapidly 
- this is the phase when the star leaves the RGB or AGB 
and moves to the left of the Hertzsprung-Russcll diagram. 
If a small companion triggers additional mass-loss while in 
the outer parts of the envelope, it will also trigger expansion 
and reduction of the envelope binding energy, because of ra- 
dius relaxation. This means that relatively light companions 
could unbind an envelope that is in principle too heavy for 
them to lift. 

The time over which stars respond to changes is the 
stellar dynamical time. For a 100-Rq RGB star this is ap- 
proximately two weeks, while for a 500-Rq AGB star it is 
of the order of a few months to a year. The latter is simi- 
lar to the CE interaction timescale determined by the few 



simulations carried o ut so far (e.g., ISandquist et al. I ll998l ; 
|Pe Marco et alj|2003l ). On these grounds we suggest that if 
the in-spiral of the companion takes longer than a dynamical 
time the giant can use its own thermal energy to help un- 
bind the envelope, but if it takes less, this energy source will 
not be available. It may be possible that when q is smaller it 
takes the companion longer to plunge into the primary enve- 
lope, thus giving the primary star more time to respond by 
expanding and reducing its own envelope binding energy. In 
other words, for slower changes the star has time to become 
virialised, which results in a stellar expansion, i.e., the star 
is tapping its own thermal energy. 

This idea is supporte d by the calculation of 
iNordhaus fc Blackmanl (|2006l . their Figure 2) that shows 
that lower mass companions take longer to in-spiral. This 
also makes sense in view of the fact that the transfer of 
orbital energy to th e envelope is via gravitational drag 
l|Ricker fc Taamll2008l ) : the companion creates a dense wake 
that slows it down and makes it fall in towards the giant's 
core. Such a wake would be less massive for lower mass com- 
panions. This suggestion can also explain the simulations' 
above-unity a values. Hydrodynamic simulations do not in- 
clude the stellar energy source, but as long as the CE in- 
spiral timescale is shorter than the thermal timescale, the 
simulated star will behave physically. 

HD 149382, with a sub-stellar companion, remains trou- 
blesome. Its value of a is so high that it cannot be explained 
simply by invoking the thermal response of the star. One 
could invoke additional sources of energy such as dissoci- 
ation and ionisation energies, but it seems suspicious that 
such source should only be invoked for the lightest compan- 
ion. It is possible that a third body was present in these 
systems, such as another planet. We could also question the 
precision of the determined parameters; however we note 
that with such a low mass companion, there is no reason- 
able way of reducing the value of a. To test this we increased 
that primary mass from 0.41 MQ(the middle of the calcu- 
lated range) to 0.47 M Q , suggested by iGeier et all (|2009l ) 
as a more likely value and we determined a under the as- 
sumption that HD 149382 is a post-RGB system as well as 
a post- AGB one. The smallest value of a (10) is obtained 
for a post- AGB object and the reduction from the value of 
33 (Table [4} is due primarily to the fact that the radius 
would be larger for a larger core mass. In our log-log fit 
this does not alter appreciably the fit parameters and the 
anti-correlation is preserved. 
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Figure 7. A plot of log(g = M2/M1) vs. log a. The simulations (crosses) are not fitted as the errors could not be determined. Top panel: 
all data is presented. Bottom panel: same as the top panel but without the post-AGB systems with no PN. The dashed lines are least 
square fits of the data. 



5 SUMMARY AND DISCUSSION 

In this paper we have compared different derivations of the 
CE efficiency equations used in the literature, and have de- 
rived a more accurate form. O ur approximation is quite sim- 
ilar to that of lWebbinkl || 1984ft . although our binding energy 
is more negative, resulting in higher values of a. 

We have then used this approximation to determine the 
value of a for a set of observed post-CE systems as well 
as a set of simulations. In so doing we have revealed the 
difficulties inherent in deriving a values from observations: 
the methods available to reconstruct the parameters of the 
primary at the time of the CE interaction do not lead to 
accurate values. As a result the errors on the determined a 
are large. 

We have also found some evidence for an anti- 
correlation between log(g = M2/M1) and log a, such that a 
«(0.05 ± 0.02)xg (_1 - 2±0 4) . Considering the exponent of q 
is < —1, this implies that smaller mass companions are left 
with comparable, if not longer post-CE periods than their 
more massive counterparts. One way this could be achieved 
is if smaller companions took longer than the primary's dy- 
namical time to penetrate into the CE than larger ones, 
giving their primary star an opportunity to react and use 
its own thermal energy to expand, thus helping to eject 
the envelope. Since the CE timescales seem to be similar 
to the stellar dynamical time scale (of the order of a year; 
ISandauist et al.lll998l ; TPe Marco et alj|2003ft this suggestion 
is plausible. 

The energetics of the CE interaction are complex, how- 
ever, and it is not enough to state that ionisation and dissoci- 
ation energies could be available, one has to also show that 
once tapped, these sources aid the envelope ejection pro - 
cess, a thing that was questioned bv lSoker fc Harpaa (|2003ft 
who argued that the opacity of the envelope would be too 
low. Several other mechanisms could play a role in the com- 
mon envelo pe ejections (e.g., excitation of pulsational waves 
ISokerl 1 1992ft . which may also depend on the in-spiral time 
scales. 



IPolitano fc Weilerl (|2007m and iDavis et all l|2010ft cal- 
culated population synthesis models under the assumption 
of a constant a, or a proportional to secondary mass. 
IPolitano fc Weilerl |2007ft compared the predicted post-CE 
period distribution, secondary and primary mass distribu- 
tions with observations of present day post-CE systems 
and cataclysmic variables. The predicted distributions are 
not strongly dependent on the prescription of a they used, 
leaving a unconstrained. Their period distribution for cur- 
rent day post-CE systems predicts many more systems at 
lo nger periods than obse rved (e.g., Miszalski et al. (2009) 
or ISchreiber et al l (|200Sf > and ISchreiber et al l i|2009ft ). The 
predicted secondary mass distribution seems also in con- 
trast with observations that show that the most represented 
companion mass has spectral type M3. 5, similarly to th e 
distribution of field main sequence stars (|Farihi et al.[ [2005). 
iPolit ano fc Weilerl (|2007ft claim that the large fraction of 
brighter companions they predict would be missed by sur- 
veys be£aus^_they would outshine the WD primary. How- 
ever [Holberg^et^al] I^OOct ) shows that the companions that 
would be missed because they outshine the primary WD are 
brighter than mid-to-late K. This means that if the real com- 
p anion mass distribution we re that predicted by the models 
of IPolitano fc Weilerl (|2007m . we would have detected more 
companions with mass larger than ~0.4 Mq. 

IDavis et all l|2010ft predict the space density of post-CE 
systems, the distribution of secondary masses and the period 
distribu tion. The lat t er tw o predictions do not match obser- 
vations. IDavis et al.l l|2010ft claim that similarly to observa- 
tions they predict a steep decline of the systems with sec- 
ondary masses <0.35 Mq. However, as we can s ee from the 
exhaustive compilation of IZorotovic et all (|2010h . the num- 
ber of systems with secondary mass belo w that threshold is 
very l arge. Finally, as is the case for the IPolitano fc Weilerl 
|2007ft models, the predicted period distribution peaks at 
longer periods and is very broad, contrary to the observed 
ones, whose distribution peaks at periods shorter than a 
day and is very narrow. The period distribution is likely to 
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become the best observation against which to test our pop- 
ulation synthesis models in the future. 

5.1 A comparison of this work to that carried out 
by Zorotovic et al. (2010) 

While this paper was in review we became aware of a publi- 
cation by Zorotovic, Gaensicke & Schreiber (2010), in which 
the authors determined the values of a for a set of post-CE 
WDs and drew conclusions from their findings in much the 
same way as we have done here. Some of their conclusions 
are similar to our own, while others differ. 

Both our and their data derive from that of 
ISchreiber fe Gansicke! < |2003h . However, they added 35 sys- 
tems that their group has newly discovered. On the other 
hand, they excluded subdwarf primaries, such as all cen- 
tral stars of PN, and the systems with low mass secondaries 
(AA Dor and HD 149382). We have insured that for all the 
obj ects in common we used the same parameters. 

IZorotovic et alj |2010l ) remarked that a values from 
simulations are in the range determined by their technique 
for the observed systems. We maintain instead that there is 
no way of comparing a values from simulations to those ob- 
tained by system reconstruction, without fully understand- 
ing what a formalism was adopted by the simulation studies 
and subject to the caveats of the simulations. This is why we 
have selected those simulations that we understand the best 
and re-determined their a values using the same formalism 
use d for the observe d syste ms. 

IZorotovic et al.l IpOlOl) compared thr ee a formalisms , 
the lYungelson et all jl993h formalism, the lWebbinkl (|l984T ) 
formalism and a hybrid form between the two formalisms 
adopted bv lHurlev et al.l (|2002r ). Differences between a val- 
ues obtained with these formalism are in li ne with our dis- 
cussio n in § 12.11 However, in contrast to IZorotovic et al.l 
(2010) we have argued that it is possible to select the best 
formalism on physical g r ound s. In the end they adopted the 
formalism of IWebbinkl | 1984), making a comparison with 
our own work strai ght forward. 

We agree with I Zorotovic et alJ (|2010l ) on the treatment 
of A, i.e., that one should calculate it according to what stel- 
lar structure one is examining. On the other hand, they did 
not disclose what values they used for those cases when they 
calculated A using stellar structure calculations, only men- 
tioning that the calculation was performed. They used the 
hydrogen abundance criterion for the core-envelope bound- 
ary, rather than the remnant envelope criterion, but that 
does not change significantly the resulting value of A. 

The primary reconst ruction technique is u sually the 
hardest part to compare. IZorotovic et alj (|2010|) used the 
stellar evolution fitting formulae of iHurlev et all < |200Cf ) to 
determine core mass, radii and luminosities as a function of 
time for a range of main sequence masses. Each measured 
primary mass (M c ) has multiple matching parameter sets 
(L, R and Mms); each set can be used to determine a. In 
this way each observed system will have a range of a val- 
ues. We, on the other hand, determined only one value of 
a, but our err or bar encompasses s imilar ranges to those 
determined by IZorotovic et al. d2010f). 

For example. IZorotovic et al.l pOlOl ) did not use the 
IFMR to determine Mms for their systems, but instead de- 
termined a range of main sequence masses which, sometime 



during their evolution, will attain the same core mass as that 
of the primary of the observed systems. This is equivalent to 
using the IFMR implied by the fitting formulae, along with 
the uncertainty in those very relations and the measurement 
errors (the error propagation is not discussed in their paper, 
nor the uncertainties in the fitting relations). In the end, 
the largest difference between t he two methods is l ikely to 
lie in the difference between the IHurlev et al] (2000) fitting 
for mulae and our MESA evolutionary sequences and fits to 
the iBertelli et all d2008l) t r acks. 

The IZorotovic et al.1 (|2010T ) reconstruction technique 
does not pre-select whether the system went through a CE 
interaction on the RGB or AGB. Indeed several of their 
systems may have gone through such interaction in either 
phase. We find that all our post-RGB systems are post-RGB 
systems also in their scheme, while those of our systems with 
M c in the range 0.47-0.55 Mq, which are post- AGB in our 
scheme (recall our criterion: M c > 0.47 Mq), have a chance 
to have gone through a CE interaction on the RGB. We 
pointed out (§ I3.4.2|) that our conclusions remain true, and 
in fact become even stronger if we eliminate the entire class 
of post- AGB stars wi th no PN, leaving the central stars of 
PN (not addressed bv lZorotovic et all d2010l)) and the post- 
RGB stars (which are post-RGB stars also in lZorotovic et akl 

After determining ranges of a values for each star, 
IZorotovic et al.l l|2010h determined the most likely values 
within those ranges by carrying out a population synthesis 
exercise. We have effectively carried out a similar exercise 
by using Baysian statistics in the determination of mass and 
radius of the primary at the time of the CE, by calculating 
the differences between Mms and Mi and between Mwd 
and M c using statistical arguments. Interestingly, the con- 
clusions drawn in the two papers are similar: a values for 
post-RGB stars are higher. Avery weak to no correlation of 
a with M2 exists (Fig. [6] middle panel). Post- AGB systems 
display a mild correlation between period a nd a, while post- 
RGB systems do not (Fig. [6] bottom panel). IZorotovic et al.l 
(2010) did not test possible correlation between q and a nor 
did they carry out a statistical treatment of their data. In 
addition, they did not list their mean a values nor the corre- 
sponding stellar parameters, making it difficult to carry out 
a more thorough comparison. 



APPENDIX A: THE CORE-ENVELOPE 
BOUNDARY AND THE VALUE OF A FOR 
DIFFERENT STELLAR MODELS AND 
EVOLUTIONARY STAGES. 

To locate the core-envelope boundary one could use a crite- 
rion based on density alone (called "Density" in Table lATj) . 
such as the radius value that minimises dm/dlogr. This cri- 
terion leads to a larger radius than any criteria based on a 
physical argument. Aside from the criterion we have adopted 
(dubbed "Remnant Envelope" in Table IAT1 see also § 12.41) . 
there are two additional criteria that lead to similar values 
of A. (i) The radius where the ab undance of hydrogen is 
0.1 (this is the criterion adopted bv lDewi fc Taurisl {2000)), 
and (ii) the largest radius where the nuclear reaction rate 
decreases below a given threshold value (just outside the 
hydrogen-burning shell; dubbed "Shell" in Tabic [All . 
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Figure Al. Stellar structures for an RGB star (left panels) which had a main sequence mass of 2M and its AGB counterpart (right 
panels; for additional model parameters see Table lAlT l . Vertical solid lines represent the core/envelope boundaries chosen according to 
the "Xjj = 0.1" criterion (leftmost vertical line), "shell" criterion and "remnant envelope" criterion (second and third vertical lines from 
the left - overlapping in the left panels) and "density" criterion (rightmost vertical line). 



Criterion: Density Xh=0.1 Shell Remnant Env. 



RGB 



M C (M Q ) 


0.314 


0.313 


0.313 


0.314 


Rc(Rq) 


0.044 


0.031 


0.032 


0.040 


X h 


0.70 


0.1 


0.23 


0.69 


£ bin /10 47 (ergs) 


-9.14 


-9.62 


-9.55 


-9.23 


A 


0.40 


0.38 


0.38 


0.39 








AGB 




M C (M ) 


0.537 


0.534 


0.534 


0.534 


i?c(R ) 


0.233 


0.037 


0.038 


0.040 


X h 


0.70 


0.1 


0.28 


0.52 


£ bin /10 47 (ergs) 


-2.05 


-3.06 


-3.00 


-2.92 


A 


0.32 


0.21 


0.22 


0.22 



Table Al. Parameters at the core-envelope boundary accord- 
ing to different criteria for the 2-Mq model on the RGB 
(Mi=2.0 M Q , R = 20 R ) and AGB (Mi=1.93 M , R = 
100 R ). 



Although the abundance criterion leads to a lower value 
of the core-envelope boundary radius, a glance at Table [ATI 
shows us that the resulting values of A are quite similar. 
There we list the values of the core mass, core radius and 
hydrogen abundance of the 2-M MESA model for two key 
evolutionary phases, when the star is half way up the RGB 
(M=2.0 M , R=20 R and r=1.0 Gyr) and at the start of 
the thermally-pulsating AGB (M=1.93 M , R=100 R and 



r=1.2 Gyr), along with the derived values of the envelope 
binding energy and of A. 
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